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INTRODUCTION 


In ECON Report 76-243-1 , we presented the mathematical formu- 
lation and the basic results of ECON's optimal decision model of 
wheat production and distribution. The present report will serve 
to document the computer programs written to implement the model. 

The programs were written in APL, an extremely compact and powerful 
language particularly well suited to this model, which makes extensive 
use of matrix manipulations. 

Chapter 1 of this report presents the algorithms used and gives 
listings of and descriptive information on the APL programs used. 

Chapter 2 gives an account of possible changes in input data. 

This report should be used in conjunction with ECON Report 76-243-1, 
which gives a complete mathematical description of the model. 


V 
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1. ALGORITHMS AND NOTATION 

1 .1 Overview and Stages of Calculations 

Essentially, the algorithms used calculate value of information as out- 
put in response to two kinds of inputs: economic parameters such as elas- 

ticities, interest rates and typical production levels; and numerical descrip- 
tions of performance of production information systems. The stages involved 
in these calculations are shown in Figure 1.1. 

For each of two production information systems, called "current system" 
and "improved system," calculations are performed in three stages, dynamic 
programming, simulation, and function evaluation. 

The first stage is dynamic programming, of which the primary output 
consists of the coefficients of the value functions. The dynamic programming 
stage requires two kinds of inputs: the economic parameters; and the statis- 

tical parameters on inventories of stored and growing crops. These statistical 
parameters are mathematically determined by the functional equations solved 
in the dynamic programming stage. Ideally, then, they would be only internal 
variables to this stage, rather than input variables. However, we have found 
no satisfactory means of solving for these variables within the dynamic pro- 
gramming calculations, so an iterative procedure is used, in which the decision 
rules determined as a by-product in the dynamic programming calculations are 
used to simulate the system (the second stage), and then the statistics collec- 
ted from the simulation are used as input for a new pass through the dynamic 
programming calculations. This alternation is continued until convergence is 
achieved. The third stage of the calculations (function evaluation) produces 
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the optimal mean present value to the United States, and the whole world, 
of the produced and distributed wheat, given a specified starting level of 
storage and growing inventories. 

After finding the optimal mean present value for both the current 
information system and the improved information system, we find the value of 
the information improvement by taking the difference of these values. 

The following (18) steps are involved in performing a complete value of 
information calculation with the model: 

1. Initialize economic parameters. 

2. Initialize parameters specific to current information system. 

3. Dynamic programming--find decision rules and quadratic and 
linear coefficients of world value functions. 

4. Simulate to find statistical parameters (current information 
system) . 

5. Check for convergence of statistical parameters. If not 
converged, update and go to (3). 

6. Find quadratic and linear coefficients of U.S. value functions. 

7. Find constant terms of value functions. 

8. Simulate to collect numerous statistics. 

9. Evaluate world and U.S. value functions for selected starting 
conditions. 

10. Initialize parameters specific to improved information system. 

11. Dynamic programming— find decision rules and quadratic and 
linear coefficients of world value functions. 

12. Simulate to find statistical parameters (improved information 
system). 

13. Check for convergence of statistical parameters. If not 
converged, update and go to (11). 
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Find quadratic and linear coefficients of U.S. value functions. 

15. Find constant terms of value functions. 

16. Simulate to collect numerous statistics. 

17. Evaluate world and U.S. value functions for selected starting 

conditions. 

18. Value of information improvement — form difference of results 

of (9) and (17). 

1 . 2 Details on Stages 

1.2.1 Initialization 

The initialization of the economic parameters consists of the following 
steps. 

Step 1 Input 

(a) m = number of periods in year. 

(b) = annual price elasticity of demand for United States 
and rest of world. This is a 2-vector. 

(c) E^ - cost elasticity of production for United States and 
rest of world by period of year. This is a matrix of shape 
m X 2. 

(d) P = average price in United States and rest of world. This 
is a 2-vector. 

(3) - average annual consumption in United States and rest 

of world. This is a 2-vector. 

(f) n = average planting in United States and rest of world by 
period of year. This is a matrix of shape m x 2. 

(g) n = number of grid points in each dimension of the state vector 
for purposes of quadratic approximation of value function. 
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(h) r = annual discount rate. 

(i) (S = average annual exports from United States to rest 
of world. 

(j) T = quadratic coefficient in cost of transportation 
function. 

(k) (jj = linear coefficient in cost of transportation function. 

(?.) D = dimension of state vector by period. This is a 

vector or length m. 

Step 2 Compute periodic demand and production function parameters, 
and discount factor. 

price - $ + 2 a X period's consumption, 
cost = 6 -I- 2 Y X period's production, 

(a) a = mP / (2CE^), a 2-vector. 

(b) 6 = P (1 - , a 2-vector. 

(c) Y.jj = P^- / (2 n.. ’ 1 = 1 , .. . , m ; j = 1 , 2 . 

Y is a matrix of shape m x 2. 

(d) 6.. = P. (1 - p-5— ) , i = 1 , ... ,m; j = 1 , 2 . 

(e) P = > discount factor for a single period. 

Step 3 Definition of period-independent part of constant arrays. 
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Step 4 Initial estimates of value function coefficients and statistical 
parameters. All are chosen by rough heuristic procedure. Values are 
refined during calculations. 

(^) ^ijk ^ i~l» ...» j~l» ...»4» k— 1, ...» 4. 

Quadratic coefficients of world value functions 

(b) -t- Q..., Quadratic coefficients of United States 
ijk 

value functions 

(c) L.. 0; i=l, ..., 6; j=l» ...» 4. Linear coefficients 

* w 

of world value functions 

(d) 0. Linear coefficients of United States value 

ij 

functions. 

55 350 0 0 

45 300 0 0 

35 40 250 250 

25 40 200 250 

15 40 150 250 

5 50 100 250 

mean values of grid point components in state space. 

10 30 0 0 

10 30 0 0 

5 1 25 5 

4 1 25 5 

3 1 25 5 

2 1 25 5 

standard deviations of grid point components in state space 
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1.2.2 Current System Initialization 

2 

Step 1 Define 2 x 10 array, Z , specifying the performance of the 
current information system. 

Step 2 If better approximations of Q, Q^, L, L^, U , and Z are 

X X 

available than those of Section 1.2.1, initialize them accordingly. 

1.2.3 Dynamic Programming 
Step 1 i = 6, period counter. 

Step 2 Define constant arrays specifying state transformation and 


incremental value functions. 

(a) New rows and columns are inserted, depending on the value of i, in 
the matrices A1 , A2, to create the three dimensional array A. The 
insertions are the period-dependent elements. The first component 
of A is used to label the two value functions (whole world and 
United States). For each fixed value j of the first component, A. 
is a symmetric matrix as follows. If i = 2 or i = 5, then 


0 


0 


0 


0 


0 
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If i = 1 , i = 3, or i = 4, then 
= Al, A2 = A2. 

If i = 6, then A^ and are just as for the case of i = 2, except the 
third row and third column are omitted. 

(b) New elements are inserted, depending on the value of i, before 

and after element 3 of B1 and B2, to create the matrix B. The first 
component of B (row index) is used to label the two value functions. 
If i = 2 or i = 5, then 

S, 6, 0 ° 1 ' 

If i = 1, i = 3, or i = 4, then B is the 2x3 matrix given by 



If i = 6, then B is just as for the case of i = 2, except that 
the third column is omitted. 

(c) The matrix C is defined, depending on the value of i, for use 
in building the quadratic programming tableau. 

If i = 1 , i = 3, or i = 4, then 

1 1 0 \ 

0 0 1 / ’ 

If i = 2, or i = 5, then 

110 0 0 
0 0 0 1 0 






If i = 6, then 



1 


1 


0 


0 


c = 


\ 0 0 


(d) State transformation matrices M and N are defined, depending on 
the value of i, as follows. 

If i = 1 , then 



\ 

/ 


\ 

1 


^ -1 

-1 



, N = 




\o 

1 / 

1 0 

0 



If i = 2, then 



/ 1 



/■’ 

-1 

0 

0 

0 


' 0 

0 ' 


' 0 

0 

1 

0 

0 

M = 

i ° 


, N = 


1 

0 

-1 

0 


\ 0 

0 1 


\ 0 

0 

0 

0 

1 

= 3 or i = 

= 4, then 









f ’ 

0 

0 



r' 

-1 

0 


f 

0 

1 

0 

0 


f 

0 

0 

0 

M = 


0 

1 

0 

1 

, N = 

1 ° 

1 

-1 


\ 0 

0 

0 

, / 


' 0 
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If i = 5, then 



r 

0 

0 

0 ^ 

1 


-1 

0 

0 

0 

0 


0 

1 

0 

0 

' 

1 

0 

0 

0 

0 

0 

0 

M = 

0 

1 

0 

1 

0 

, N = 

i “ 

1 

0 

0 

-1 

0 


^0 
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0 

1 i 

f 

\ 0 

0 

0 

0 

0 

1 

If i = 6, 

then 
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step 3 Build a matrix X whose columns are the grid points for 
value function approximation. X has D^. rows and columns. In 
each coordinate j of the state space (of dimension ) , a set 'Sj of 
n equally spaced points in formed with mean U . and standard devi- 

X j j 

ation E .. Each column of X is formed from one of the combi- 

X J J 

nations of elements, the first element from , the second from 
, . . . , the D . th from . 

2 , 1 U*! 

step 4 j = 1, counter for state points. 

Step 5 Quadratic Programming 

(a) S = column j of X, the state point. 

(b) E = + p N' N. 

(c) F = + 2p MSN + p N. 

(d) G = p(S'M' MS + MS). 

(e) ( S if D. = 2 
D = I 

( (s^, s^) if D. = 4 

(Note on notation. When i = 6, we take i + 1 reduced mod 6. 

That is, for this case, i + 1 = 1. When = 2, the matrix 

Q-j+] is 4 X 4, but only the upper left block of size 2x2 

is significant--the other elements are never used nor changed 

during the algorithm. For this case, we write "Q^+i" to 
mean only this 2x2 matrix. Similarly, means only 

the first elements of the 4-vector.) 

(f) Find 1-stage decision vector Y to maximize Y'E Y + Y’F 
subject to the constraints 
Y > 0 , CY < D . 
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step 6 Evaluate the cumulative value function for period i at 
point S. 


V, = VE ,,, Y + 

Step 7 End of loop on state points. 

(a) j ^ j + 1 

(b) Go to Step 5 if j < , Step 8 otherwise. 

Step 8 pind period-i coefficients of cumulative value function 
by least squares approximation. 

(a) Form the matrix MAT which, when postmultiplied by the 
values of a function at the points which are the columns 
of X, produces the coefficients of the least squares fit 
to the function at those points. 


(b) COEFF = (MAT)V 

(c) Select Q and L from COEFF 
Step 9 End of loop on periods. 

(a) i i - 1 

(b) Go to Step 2 if i > 1 , Step 10 otherwise. 

Step 10 End of loop for convergence. 

Go to Step 1 if Q and L are unchanged from last iteration; otherwise 
stop. 


1.2.4 Simulation 

The simulation algorithm is an adaptation of the dynamic programming 
algorithm. It is run after the value function coefficients Q and L have 
been calculated, and uses many of the same variables. 

Step 1 Initialization 

(a) Input Z array, a 2 x io matrix giving the means of the 

squares of the stochastic term <J>^- of the state transformation. 
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(b) Input n, the number of years of simulation to be performed. 

(c) Create standard deviation array SD for use in sampling 
values of <I>^. by reshaping the Z array and taking square root. 

In its original shape, the first row refers to the Untied States, 
and the second to the rest of the world. Each row contains terms 
referring to the growing crop before the beginning of the crop 
year, followed by terms referring to the same crop during the 
crop year, considered to begin June 1. The reshaped array has 
shape 4 X m; the growing crop and current crop at a given time 
of year are represented by different rows in the same column. 
Formally, the process is as follows. 


1 1 1 

1 1 2 

2 2 0 

2 2 0 / 

9 10 5 

3 4 5 

9 10 0 

3 4 0 I 


SO. . = 


where 


T = 




1 

2 

0 

0 


1 

1 

2 

2 


1 

1 

2 

2 


U = 


6 

0 


8 

2 

8 


0 1 


2 

and Z is taken to be 0. 

0,0 
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(d) Select vector $ with 4nm elements by random sampling. For each 
of the mn periods of the simulation, four successive elements of 

4> are obtained from discrete uniform distributions with mean 0 and 
standard deviations SD, i, where i is the index (1 to 6) of the 

current period. 

(e) Find standard deviations SD' over 2n years of the blocks of 4m 
adjacent elements of the vector (3>, - $). 

(f) SD ^ SD^ / SD' 

(g) Repeat step (d). At this point, the sample points comprising 
(<I), - <I>) have mean 0 and the standard deviations of the original 
SD array. 

(h) S = first two components of first row of U , the mean inventory 

A 

level vector at time 1. 


(i) 


U 


X 


z 

X 


/ 0 0 0 
0 0 0 
0 0 0 
0 0 0 
10 0 0 
1 0 0 0 


0 
0 
0 
0 
0 
0 / 


initial ize 


mean and standard deviation of state vector. 

Step 2 t ■<- 0, counter fok' vector. 

Step 3 i 6, period counter. 

Step 4 Increment grid parameter accumulators. 

C 

(a) U . U ■ + zero components to S). 

X5I Xjl fcil 1 

S2 

^x,i ^x,i 2n ^i ” components to S). 
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Step 5 Define constant arrays specifying state transformation and 
incremental value functions. This step is identical to Step 2 of 
Section 1.2.3. 

Step 6 Quadratic Programming. This step is identical to Step 5 of 
Section 1.2.3, except that 5a is omitted, since S has already been 
selected. 

Step 7 Apply state transformation to get next period's state vector. 

(a) S 'f- MS + NY + ‘*’£+3’ 

(b) fi, ^ J, + 4 

Step 8 End of loop on periods. 

(a) i f- i + 1 

(b) Go to Step 4 if i < 6, Step 9 otherwise. 

Step 9 End of loop on years. 

(a) k ^ k + 1 

(b) Go to Step 3 if k < n. Step 10 otherwise. 

Step 10 End of loop on antithetic variates. 

(a) 4) - 4) 

(b) Go to Step 2 if Jl < 8 mm. Step 11 otherwise. 


Step 11 Complete statistical calculations. 


n(r, - (U^)^) 

n - 0.5 


Step 12 Stop. 
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1.2.5 United States Value Functions 

The algorithm of this section is an extension of the dynamic pro- 
gramming algorithm (Section 1.2.3) to calculate the United States value 
function, as well as the world value function which is maximized. 

Steps 1 through 5 Identical to Steps 1 through 5 of Section 1.2.3. 
Step 6 Evaluate the cumulative value functions (world and United 
States) for period i at point S. 

Vi - 

= A2 + pN'Q'I^^N . 

f"" = B2 + 2pq‘J^^msn + plV^^n . 

g"' = p(S'M'Q'j'^^MS + lIJ^^MS) . 

yU ^ y,^u Y ^ y'fV + gV 
1 1+1 1+1 1+1 

Step 7 End of loops on state points. 

(a) j ^ j + 1 

(b) Go to Step 5 if j < n^'' , Step 8 otherwise. 

Step 8 Find period-i coefficients of cumulative value function 
(world and United States) by least squares approximation. 

(a) New rows and columns are inserted, depending on the value 
of i , in the matrices A1 , A2, to create the three dimen- 
sional array A. The insertions are the period-dependent 
elements. The first component of A is used to label the 
two value functions (whole world and United States). For 
each fixed value j of the first component, A. is a symetric 

s) 

matrix as follows. If i = 2 or i = 5, then 
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If i = 1 , i = 3, or i = 4, then 
= A1 , = A2. 

If i =6, then A^ and A^ are just as for the case of i =2, 
except the third row and third column are omitted. 

(b) COEFF ^ (MAT)V 

(c) COEFFU ^ (MAT)V^ 

(d) Select Q, L, and KEF from COEFF; L^, and KEFU from 
COEFFU. 

Step 9 End of loop on periods. 

(a) i i - 1 

(b) Go to Step 2 if i >1, Step 10 otherwise. 

Step 10 End of loop for convergence. Go to Step 1 if any of Q, L, 
Q^, are different from last iteration; otherwise stop. 
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1,2.6 Constant Terms 

The constant terms of the cumulative value functions {world and 
U.S.) are most efficiently calculated in this separate algorithm rather 
than in the dynamic programming algorithm since convergence is slow and 
the results are not needed for the determination of the decision rules, 
but only for evaluation of the cumulative value functions. 

Step 1 Initialization. Identical to Step 1 of Section 1.2.4., 
except (b) is omitted. 

Step 2 i = 6, period counter. 

Step 3 Update estimate of period-i constant terms of cumulative 
value functions (world and U.S.), 




M K, . KEF. . p * E 

mm sj ^ • 


(b) kV ^ KEFU. + p 




(Note on notation. When i = 6, we take i + 1 reduced mod 6, as 
explained in the note after Step 6 of Section 1.2.3.). 

Step 4 End of loop on periods 

(a) i i - 1 

(b) Go to Step 3 of i >1, Step 5 otherwise. 

Step 5 End of loops for convergence. If K or has changed since 
last iteration, go to Step 2; otherwise Step 6. 

Step 6 Stop. 
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1.2.7 Statistics 

The statistical summaries of the operation of the wheat markets are 
obtained by simulation; this alogrithm is thus essentially the same as 
that of Section 1.2.4, but additional quantities are tracked. 

Step 1 Initialization 

(a) Input Z array, a 2x10 matrix giving the means of the 
squares of the stochastic terms of the state transfor- 
mation. 

(b) Input n, the number of years of simulation to be performed. 

(c) Create standard deviation array SD for use in sampling 
values of by reshaping the Z array and taking square 
root. In its original shape, the first row refers to the 
United States, and the second to the rest of the world. 

Each row contains terms referring to the growing crop 
before the beginning of the crop year, followed by terms 
referring to the same crop during the crop year, considered 
to begin June 1. The reshaped array has shape 4xm; the 
growing crop and current crop at a given time of year are 
represented by different rows in the same column. Formally, 
the process is as follows. 




^ij’^'ij 
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where 



and Y? „ is taken to be 0, 

0,0 

4 Columns 9 Columns 12 Columns 



multiplier to weight tracked quantities for annual mean, 
total, or present value calculations. 

(e) Select vector $ with 4nm elements by random sampling. For 
each of the mn periods of the simulation, four successive 
elements of <J) are obtained from discrete uniform distribu- 
tions with mean 0 and standard deviations SD ., where i is 
the index (1 to 6) of the current period . 
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(f) Find standard deviations SD' over 2n years of the blocks 
of 4m adjacent elements of the vector (4>, -4>). 

(g) SD ^SD^/SD' 

(h) Repeat Step (d). At this point, the sample points com- 
prising ($, -$) have mean 0 and the standard deviations 
of the original SD array. 

(i) S = first two components of first row of U^, the mean 
inventory level vector at time 1. 


(j) 


^0 0 
0 0 


y = VAR = 


0 0 
0 0 
0 0 


0 0 



(27 columns). 


initialize mean and variance by period of 27 quantities to 
be tracked for statistics. These are: 

(1) U.S. Stock Estimate 

(2) U.S. Growing Crop Estimate 

(3) R.O.W, Stock Estimate 

(4) R.O.W. Growing Crop Estimate 

(5) U.S. Consumption 

(6) U.S. Exports 

(7) U.S. Planting 

(8) R.O.W. Consumption 
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(9) R.O.W. Planting 

(10) Random Change in U.S. Stock Estimate 

(11) Random Change in U.S. Growing Crop Estimate 

(12) Random Change in R.O.W. Stock Estimate 

(13) Random Change in R.O.W. Growing Crop Estimate 

(14) U.S. Price 

(15) R.O.W. Price 

(16) U.S. Export Revenue 

(17) U.S. Production Cost 

(18) R.O.W. Production Cost 

(19) R.O.W. Transportation Cost 

(20) U.S. Gross Welfare 

(21) R.O.W. Gross Welfare 

(22) World Net Welfare (all agents) 

(23) U.S. Net Welfare (all agents) 

(24) U.S. Consumers' Net Welfare 

(25) U.S. Producers' Net Welfare 

(26) R.O.W. Consumers' Net Welfare 

(27) R.O.W. Producers' Net Welfare 

(k) AVAR (0, 0, ..., 0), initialize annual means and 


variances of above 27 quantities. 

(1) Input S, the state vector at the start of the simulation, 
time 1 . 

Step 2 £ ■*- 0, Counter for $ vector. 


Step 3 i ■<- 6, Period counter. 

Step 4 Define constant arrays specifying state transformation and 
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incremental value functions. This step is identical to Step 2 of 
Section 1.2.3. 

Step 5 Quadratic Programming. This step is identical to Step 5 of 
Section 1.2.3, except that 5a is omitted, since S has already been 
selected. 

Step 6 Build vector STAT^ of 27 quantities at time i for statistics. 

(a) STAT. , ..., STAT.. are taken from S. 

(b) STAT^.^, ..., STAT^g are taken from Y. 

(c) (STAT^.^q, .... STAT^.^^) •••, <i>£+ 4 )- 

(d) (STAT.^^, STAT.^ 3 ) - ( 6 ^ + 2a^ STAT.^, 6 ^ + 2a^ STAT.g), 
prices. 

(e) STAT. STAT.,^ x STAT._ (export revenue). 

(f) (STAT.^^, STAT.^g) - (STAT.^, STAT.^) x 
(6. + (STAT, STAT^g)Y^-)» production cost 

(g) ^TAT^ STAT^g x (w + STAT^.^ x x), transportation cost. 

(h) (STAT.^g, STAT- 2 ^) - (STAT .5 X (B^ + x STAT.^), 

STATis X (^2 “2 9 >^oss welfare 

(i) (STAT^. 22 > •••. STAT.^.^) ^ Y'AY + Y'B, net welfare. 

Step 7 Increment Statistical Accumulators. 

STAT. 

^ ^ * TT ■ 

(STAT 

VAR. ♦ VAR. + — 5 -!- . 

Step 8 Apply state transformation to get next period's state vector. 

(a) S - MS + NY + (<j)^^^, ..., 4 .^^^). 

(b) £ -<- + 4 . 
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Step 9 

(a) 

(b) 
Step 10 

(a) 


(b) 


Step 11 

(a) 

(b) 
Step 12 

(a) 

(b) 
Step 13 

(a) 

(b) 


Step 14 
Step 15 


End of loop on periods, 
i i + 1 

Go to Step 4 if i < 6, Step 10 otherwise. 


6 

V STAT. X mult. 

1 1 


AVAR ^ AVAR + 


(E STAT. X mult.) ^ 
i=l ^ ^ 


calculate annual means and mean squares of tracked 
quantities. 

End of loop on years, 
k ^ k + 1 

Go to Step 3 if k < n. Step 12 otherwise. 

End of loop on antithetic variates. 


4 > ■*- - $ . 


Go to Step 2 if X, < 8mn, Step 13 otherwise. 


VAR ^ ~ , 

(AVAR = vhn 

AVAR ^ i— ^ , 

n — 

complete calculation of periodic and annual variances of 
tracked quantities. 

Print u, VAR, and AVAR arrays. 

Stop. 
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1.2.8 Evaluations 

This algorithm evaluates the quadratic functions giving optimal mean 
net present value to the whole world and to the United States, for chosen 
arguments (starting inventory estimates). As an option, the algorithm 
calculates the mean value of these value functions over the various pos- 
sible starting inventory arguments. 

Step 1 Input argument S = (s^, s^) 

Step 2 If S = 0, go to Step 5. 

Step 3 


(a) Print r[K^ + S'(L^ + Q^^S)], whole world value. 

(b) Print r[K^ + S'(L^ + Q^S)], United States value. 
Step 4 Go to Step 6. 

Step 5 


(a) Print t S'(L, t Q^s)]. 

whole world mean value. 

(b) . Print * kJ . S'dlJ . Q^S)]. 

United States mean value. 

Step 6 Stop. 

1.2.9 Value of Information 

The annual value of the improvement in information in the steady 
state is obtained as follows. 

Step 1 Do the evaluation of Step 1.2.8 (mean value case) with K, 

L, Q, K^, L^,andQ^ corresponding to the current information system 
and S obtained from a simulation corresponding to the current infor- 
mation system (S is defined in Section 1.2.4, Step Ih). Denote the 
results (W^, V^). 
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Step 2 Do the evaluations of Step 1.2.8 with K, L, Q, K^, L^, and 
corresponding to the improved information system and S obtained 
from a simulation corresponding to the improved information system. 
Denote the results (Wj, Vj). 

Step 3 The value of the improvement to the world is Wj - W^, and 
to the United States in Vj - V^, 

Step 4 Stop. 

1 . 3 The APL Programs 

The dynamic programming algorithm, which is the heart of the model, 
is programmed in the APL function VFS (standing for value functions). 

The initializations required are handled by the function INITIAL. The 
simulation to find the statistical parameters on the state variables is 
done by SIMGRID, and the alternation between VFS and SIMGRID required to 
obtain convergence of the statistical parameters is performed by ITERATE. 
Thus, the user begins a calculation by running INITIAL and ITERATE, with 
VFS and SIMGRID called automatically from ITERATE. The quadratic and 
linear coefficients of the United States value functions are obtained 
•after convergence by USVFS, and the constant terms of the value functions 
are calculated by CONTERMS. The statistics are obtained by simulation 
in SIMSTATS, and finally, the evaluation of the quadratic value functions 
is done by EVAL. 

Besides the functions mentioned above, there are various subordinate 
functions which they call. The following sections contain a glossary of 
the main global variables, and a listing and discussion of each APL 


function. 



1,3.1 Global Variables 


With each APL variable, we give the shape, algebraic name used in 
this report, a brief description, and the step in the algorithm where 
the variable is defined. 


APL Name 

and Shape 

A1 qebraic 

Name and Description 

A, 2 4 4 


A 

Quadratic coefficients of incre- 
mental value functions. Section 
1,2.3, Step 2a. 

ADR, 1 


r 

Annual discount rate. Section 
1.2.1, Step Ih. 

A, 2 4 4 


A 

Quadratic coefficients of incre- 
mental value functions. Section 
1.2.3, Step 2a. 

AFLAG, 1 


— 

Flag used in function NOTESHIFT. 

AIN, 2 3 

3 

A1 , A2 

Period-independent part of A, 
Section 1.2.1, Step 3. 

ALPHA, 2 


a 

-h Slope of demand function. 
Section 1.2.1 , Step 2a. 

APRD, 2 


— 

Average annual production. 

B, 2 4 


B 

Linear coefficients of incremental 
value functions. Section 1.2.3, 
Step 2b. 

BETA, 2 


6 

Intercept of demand function. 
Section 1.2.1, Step 2b. 

BIN, 2 3 


Bl, B2 

Period-independent part of B, 
Section 1.2.1, Step 3. 

C, 2 4 


C 

Matrix of left hand side coef- 
ficients in quadratic program- 
ming problem. Section 1.2.3, 
Step 2c. 

CNSMPTN, 

2 


Average annual consumption. 
Section 1.2.1, Step le. 

CUR, 2 10 


Description of current informa- 
tion system. Section 1.2.2, 

Step 1 . 



APL Name and Shape 


Algebraic Name and Description 


DELTA, 6 2 

6 

Intercept of production cost 
function. Section 1.2.1, Step 2d 

DIMX, 6 

D 

Dimension of state vector by 
period. Section 1.2.1, Step U. 

UJ 

E 

Quadratic coefficients in expan- 
sion of value functions in powers 
of Y. Section 1.2.3, Step 6b. 

EXPRTS, 1 

e. 

Average annual exports from United 
States to rest of world. Section 
1.2.1, Step le. 

F, 4 

F 

Linear coefficients in expansion 
of value functions in powers of 
Y. Section 1.2.3, Step 6c. 

Ganma , 6 2 

Y 

Slope of cost function for pro- 
duction. Section 1.2.1, Step 2c. 

GE, 2 10 


Description of improved information 
system. 

IMAX, 1 

• • 

Number of state points at which 
value functions are evaluated, 
n^i . Section 1.2.2, Step 8b. 

KEF, 6 

KEF 

Part of the constant terms of the 
world cumulative value functions. 
Section 1.2.5, Step 9c. 

KEFU, 6 

KEFU 

Part of the constant terms of the 
U.S. cumulative value functions. 
Section 1 .2.5, Step 9c. 

KON, 6 

K 

Constant terms of the world cumula- 
tive value functions. Section 1.2.6, 
Step 3a. 

KONU, 6 

K^ 

Constant terms of the U.S. cumula- 
tive value functions. Section 1.2.6, 
Step 3b. 

L, 6 4 

L 

Linear coefficients of world cumula- 
tive value functions. Section 1.2.3, 
Step 9b, and Section 1.2.1, Step 4c. 

LU, 6 4 


Linear coefficients of U.S. cumula- 
tive value functions. Section 1.2.5, 
Step 9c, and Section 1.2.1, Step 4d. 



APL Name and Shape Algebraic Name and Description 


M, 2 2 
4 2 

M 

State transformation matrix. 
Section 1.2.3, Step 2d. 

4 4 
2 4 

MEANX, 6 4 

U 

X 

Mean values of grid point com- 
ponents in state space. Section 
1.2.1, Step 4e. 

N, 2 4 
2 3 

N 

State transformation matrix. 
Section 1.2.3, Step 2d. 

4 5 
4 3 
4 6 

NGRID, 1 

n 

Number of gird points in each dimen- 
sion of the state vector. Section 
1.2.1, Step Ig. 

NPERIODS, 1 

m 

Number of periods in year. 

OMEGA, 1 

w 

Linear coefficient in cost of trans- 
portation function. Section 1.2.1, 
Step Ik. 

PRDCTN, 6 2 

n 

Average planting in United States 
and rest of world by period of year. 
Section 1.2.1, Step If. 

Q, 6 4 4 

Q 

Quadratic coefficients of world 
cumulative value functions. Section 
1.2.3, Step 9b., and Section 1.2.1, 
Step 4a. 

QU, 6 4 4 


Quadratic coefficients of U.S. 
cumulative value functions. Section 
1.2.5, Step 9c, and Section 1.2.1, 
Step 4b. 

QPREV, 644 

-- 

Stored value of Q for convergence 
check. Section 1.2.3, Step 11. 

QPREVU, 644 

— 

Stored value of Q^ for convergence 
check. Section 1.2.5, Step 11. 

RHO, 1 

p 

Discount factor for one period. 
Section 1.2.1. Step 2e. 

STDX, 6 4 

E 

X 

Standard deviations of grid point 
components in state space. Section 
1.2.1, Step 4f . 



API Name and Shape Algebraic Name and Description 
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TAU, 1 

T 

Quadratic coefficient in cost of 
transportation function. Section 
1.2.1, Step 1 j . 

TOL, 1 

— — 

Tolerance factor for convergence 
of value function coefficients. 
Defined from keyboard. 

X, 2 9 
4 81 

X 

Matrix of grid point coordinates. 
Section 1.2.3, Step 3. 

Y, 3 
4 

Y 

1 stage decision vector. Section 
1.2.3, Step 6f. 


5 

1.3.2 Function INITIAL 

Function INITIAL is listed in Figure 1.2. If a new calculation is 
desired in which demand or cost parameters are changed, or the discount 
rate, or the degree of resolution (NGRID), the appropriate changes 
should be made in lines [2] through [13] of INITIAL. Then INITIAL 
should be run, followed by VFS or ITERATE. 

INITIAL is quite straightforward, but some explanation of line [6] 
is in order. Planting of wheat is assumed to take place in the United 
States in the second and fifth periods, and in the rest of the world 
in the second, fifth, and sixth periods. Thus, the production cost 
function parameters, y and 6 , are used only in these periods (by 
function STRUCTURE). The elements of the array PRDCTN corresponding 
to the non-planting times, therefore, have no effect on the calculation, 
but if they are zero, numerical problems are encountered in line [16]. 

To avoid this, we have put 0.1 in these positions. 

1.3.3 Function ITERATE 

This function, listed in Figure 1.3 performs the iterations of dynamic 
programming and simulation leading to convergence of the statistical 
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Figure 1.2 Listing of APL Function INITIAL 
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Figure 1.3 Listing of APL Functions ITERATE and NEWGRID 
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2 

parameters. Its left argument is the array Z (Section 1.2.4, Step la) 
describing the performance of the information system. Its right argument 
is 1/2 the number of years to run the simulations. After running SIMGRID 
to determine the simulated statistical parameters, ITERATE calls NEWGRID 
to form the new parameters as a linear combination of the old ones and 
the simulated ones. The relaxation coefficient, STEP, is a global vari- 
able defined from the keyboard. Values of STEP between .05 and .5 have 
been found suitable. The variable NUM is also defined from the keyboard 
and is used as the argument of VFS to limit the number of dynamic pro- 
gramming cycles before SIMGRID is rerun. NUM=1 has been found most 
often suitable for speedy convergence. 

1.3.4 Function VFS 

The listing of this function is in Figure 1.4. The following is a 
glossary of important variables used in VALUEFUNCTIONS and not discussed 
in Section 1.3.1. 

APL Name and Shape Algebraic Name and Description 


CLOOP, 1 


Label for start of loop for 
convergence of value function 
coefficients. 


CDEFF COEFF 

(Vector, length varies) 


Receives results of least squares 
approximation of value function 
coefficients. Defined in Function LS. 


DIM, 1 


Right hand side of constraint 
inequality on decision vector Y 
in quadratic programming problem. 

Dimension of state variable in 
next period. 

Constant term in expansion of 
value functions in powers of Y. 

Counter for state points. 
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Figure 1,4 Listing of APL Function VFS 



APL Name and Shape 

Algebraic 

Name and Description 

INEXT, 1 

i+1 

Index of next period of the year. 

IPER, 1 

i 

Counter for period of the year. 

ITERS, 1 

-- 

Maximum number of iterations in 
seeking convergence. 

K, 1 

— 

Counter for iterations in seeking 
convergence. 

LLTD, 4 
2 

L 

Row of L referring to next period 

NC, 1 

— 

Number of columns in quadratic 
programming tableau. 

NCY, 1 

— — 

Dimension of decision vector, 
number of columns in constraint 
matrix C. 

NR, 1 

— 

Number of rows in quadratic 
programming tableau. 

NRY, 1 

— 

Number of rows in constraint 
matrix C. 

PERLOOP, 1 

-- 

Label for start of loop over 
periods of year. 

QLTD, 4 4 
2 2 

Q 

Submatrix of Q referring to next 
period. 

R 

(vector, size varies) 

— 

Indices of basic columns in qua- 
dratic programming algorithm. 

S, 2 
4 

s 

State point. 

STATELOOP, 1 

— 

Label for start of loop on state 
points. 

VAL, 9 
81 

V 

Vector of values of cumulative value 
function at state points. 


The function VFS is a direct implementation of the dynamic program- 
ming algorithm, as presented in Section 1.2.3. It calls the subordinate 
functions DRAW, STRUCTURE, BUILDX, TABLEAU, MAXIMIZE, LS, and CONVERGE. 



DRAW is called at the beginning of execution to select the random 
sample of values of . 

STRUCTURE is called at the beginning of the loop on periods of the 
year to set up the arrays A, B, C, M and N, Thus, it covers Step 2 of 
Section 1.2.3. 

BUILDX covers Step 3, the formation of the array of grid points. 

TABLEAU prepares the initial tableau for the quadratic programming 
calculations, which are actually carried out in MAXIMIZE. 

LS performs the least squares fit to obtain new Q and L. 

CONVERGE compares the new Q array with the one stored in the last 
iteration, and takes the value 1 when convergence has been achieved. 

1.3.5 Function SIMGRID 

Figure 1.5 gives a listing of this function. The following is a 
glossary of important variables used in SIMGRID and not discussed above. 


APL Name and Shape 

Algebraic 

Name and Description 

ANTLOOP, 1 

— 

Label for start of loop on 
antithetic variates. 

FILLP, 6 4 

-- 

Array of flags to expand 4> to 
maximum dimension of state space. 

FILLX, 6 4 

-- 

Array of flags to expand S to 
maximum dimension of state space. 

KNT, 1 

“ — 

Pointer for array of random terms. 
Locates starting point for this year 
and period. 

ND, 1 

-- 

Dimension of state space next 
period, number of rows of M. 

PH (vector, length 

— 

Random terms selected in DRAW. 


depends on ITERS) 


PHI, 2 
4 


Random term in state transformation. 
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v!r>xMC-ir:;ir'[;0]'? 

V SXGir>« SIMGF.;TI> ITERG 

Ci: 6 4 f 1 0 1 0 1 0 1 0 vl6fl 

C2: r.Ti...LR<.. 6 4 f( 10 10 )f(16fl)y 10 10 

C 33 1 >'■ \ (. iters , . S.TGSQ ) rz ( STTERS ,, , SSITtSC ) 

[;4;] DRAW 

C 5 3 I- 1 : 

r 6 3 T s r A T E i- 2 1 M E a r - j ; ; |;; j ; "j 

C7;] MEAM?;<-GTi:iKf. <fj 4 f-r. 

He 3 ?3GHz-.l 

i;; 9 ;;| am t !„. O O *. K T f ■ (> 

r 1 0 3 - ^ f ' 

ril 3 S^-STSTATE 

ri 23 i;fi 

ri33 tearloop j 

[;143 I ERR El 

Cl 5:1 PERI... OOP: 

C 1 6 3 ^ [■ I E R 5 "I <- M E A N [; x f - e n j -j .}. ( 4 f s ) 1 t e p: s x 2 
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C 25 3 
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Figure 1.5 Listing of APL Function SIMGRID 



SGN, 1 


Multiplier (+1) to fix sign of 
random terms to achieve antithetic 
variates. 

SIGSQ, 2 10 Input array describing performance 

of information system. See 
Section 1.3, Step la . 

XB, 4 X State variable expanded to maximum 

dimension . 

The function SIMGRID is a direct implementation of the simulation 
algorithm presented in Section 1.2.4. It uses some of the same subor- 
dinate functions as VFS, namely STRUCTURE, TABLEAU, and MAXIMIZE. It 
uses the subordinate function DRAW to select the random terms. 

1.3.6 Function USVFS 

The listing of this function is in Figure 1.6. It does all the 
calculations of VFS, and simultaneously calculates QU, LU, KEF, and 
KEFU, as discussed in Section 1,2.5 For the least squares fit, USVFS 
calls the function LSU, and for the convergence check, it calls 
CONVERGEU; otherwise, it uses the same subordinate functions as VFS. 

1.3.7 Function CONTERMS 

This function is listed in Figure 1.7. It is a straightforward 
implementation of the algorithm described in Section 1.2.6. 

1.3.8 Function EVAL 

This function is listed in Figure 1.7. It is a straightforward 
implementation of the algorithm described in Section 1.2.8. 

1.3.9 Function SIMSTATS 

Figure 1.8 gives a listing of this function. The following is a 
glossary of important variables used in SIMSTATS and not discussed 
above. 
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VUSVFSCD]^ 

V SIGSa USVFS ITERS 

Cl] FiLLKf. 64fl0101010 rl6Fl 

123 FILLF<- 64fl01011111111111111111010 

C33 K^l 
CA 3 cloof: 

[5] IFER «-6 
C63 ferloof; 

C7] structure: ifer 

eg] IFER FUIUI'K WGRIII 

C93 IFI 

C lOD IWE::Tf.l+WF ERIOnSxIF ER 

Cll] X) [23 

[123 t'lMf r>IMX[ IXEXT] 

[133 auTfif. < r<iM , r<iM ) f «[ ihext j j ] 

[14] L.l.TD4-EiIMf L[ IMEXT ; ] 

C153 E:<-A[1 J f ] + RHOx (6)K) + . XOL.TD4-, xn 

C163 QLTnUf- ( IrlM , i;>IM ) f QU[ IMEXT J J ] 

[17] LLTr'U«-r.IMfLU[If<EXT; ] 

Cl 83 e:u<-aC2» ; ]+f;hox ( JsiN) + ♦ x«LTr>u+ * x« 

[19] MCT<-(fC)[23 

C203 »<P'»XfC)[i] 

C213 TABLEAU 
[223 HP:F(fTAB)[l] 

C233 HC>l+(f TAB) [23 

[24 3 Rf-HC t + \ HF 

[253 VALflMAXfO 

C263 VALU<-VAL 

[273 STATELOOFJ 

C283 Sf.BIMX[IFER]tX[ J I] 

C293 F"fE<[1 J 3+FHOx ( ( 2xoltb+ , x (ms«-m+, xs ) )+lutb> 4 , xx 

C303 04-RHOX (UUTri 4 -OLTr. 4 - ♦ XMS) 4 - . XMS 
C313 fufb[ 2 J ] + rhox ( ( 2xaLTriU4- , xms ) 4-l.lti;.u ) 4 » xx 
C 323 GU 4 .RHOX < LLTBU 4 -CU.Ti:tU 4 - , XMS ) 4 - ♦ XMS 

C333 xb4.fil.lx[ifer; ]\S 

[34] i.4XB[i 3] 

C353 r-FO. 00001 rr- 

C363 MAXIMISE 

[37 3 VAL[ I]«-G4- (F4-E4-, X T ) 4-, XT 
[38 3 VALU[I]f.GU 4 (FU 4 EU 4 , X'T) + , X T 

[393 - 4 STATEL 00 FX ) imax_>_im 4-1 
[40 3 EiIM<-IHMX[ IF ER] 

[41] LSUS 

C423 u<-( ,Vo . <.Vf.\niM) 

C433 T<-(BIM,BIM)pU\<-DIM4-l) 4 ,COEF-F- 

[44] Q[Ifer; (i(fT)[i]);o ( f t ) [23 ) 3<0 . 5x T4 ts)T 

C453 l[ifer; ) f T] 4 .T<-tnMt < -bim+i ) f coeff 

[46] KEF[IFER]4-ltCOEFF 

[47] TUF<BIM,rHM)fU\(-i:iIMfl) 4 ,COEFFU 

[483 QU[IFEr; ( i (fTU)[i]> ; ( \ < p tu ) [2 3 ) 3 <0 . 5x tu 4 6 jtu 

[49] LU[ IFER ; ) pTU]FTUFBIMt <-BIM4-1 ) fCOEFFU 

C503 KEFU[ IFER]F“lf COETFFU 

C51 3 - 4 FERLOOFX ) 1 <.IFEP:f.lFER~l 

C523 -»( (COHVERGEU K)=l)/0 

C533 -»CUOOFx \ ITERS 2 .K 4 .K 4 1 


Figure 1.6 Listing of APL Function USVFS 



VCONTERMSCO!]^ 

V SIGSO CONTEP:M5 ITERS 

Cl 3 Sr<2<-LOSS5H«F E RESHAPE SIGS« 

C23 »<<-i 

C33 CLOOF-J 
t;4 3 iPER<-6 
C53 REREOOPJ 

[63 tiiM<-DiM>:[ imekt«-i+6tIF e:f:] 

[7] tiSD<-riIM^5t'2C ; IF ER] 

C83 Uf.,(iDiM)o.s=ir.iM 

[93 aLTt>«- ( I'lM , DIM ) f G[ IHEXT f J 3 
[103 QLTDU<- ( DIM , DIM ) f OU[ IHEXTJ } 3 

[113 KOH[ IRER3f KEF[ IPER3 + RHOX l<OK'[ I MEX T 3 +DSD + 

C 123 KONU[ IPER3<- KEF-U[ IPER3+RHOX KO(fU[ IK'EXT3 + Dii 
C 133 -»PERLOOPx I 1 <.IPERf-IPER-l 
C143 -»CLOOPx \ ITERS2.K4-K + 1 
C153 12 2 tAI'Pxkoh 
C16 3 12 2 tAI-'Rxkohu 

V 


veval[03v 

V EVOL ORG 

[13 adrxkon[13+arg+. XL[ 1 J3+«PG+, xo[l ; ?3 
[23 ODRxKOI<U[13+«RG+. Xl-U[l J 3 + ARG+, x«u[i ; J 3 

V 


, xU/,OLTD 

♦ XU/,QLTDU 


Figure 1.7 Listings of APL Functions CONTERMS and EVAL 



VSIMSTATSi; ;rj]V 

<7 SIGSa SIMSTATS ITEP:S 

C13 MULTf.fi} (/| 6 f«VG),|:i](9 6 fl)»Ci:J(2 6 F avg<-i~ 6 ) » C 1 3 12 6 fP:HO*}6 

C23 riLDCf. 6 4 F 1 0 1 0 1 0 1 0 rlAFl 

C3] FiLUif. 6 5F1101011111 1101011010 11111 
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Figure 1.8 Listings of APL Function SIMSTATS 
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Figure 1.8 Listings of APL Function SIMSTATS (continued) 



APL Name and Shape Algebraic Name and Description 


AMN, 10 


^A 

Annual means of 27 quantities 
listed in Section 1.2.7, Step Ij. 

AVAR, 10 


AVAR 

Annual variance of 10 quantities 
listed in Section 1.2.7, Step Ij. 

FILLP, 6 

4 

— 

Array of flags to expand 4 to 
maximum dimension of state space. 

FILLX, 6 

4 

-- 

Array of flags to expand S to 
maximum dimension of state space. 

FILLY, 6 

5 

— 

Array of flags to expand Y to 
maximum dimension of decision space 

GRS, 2 


STAT(i20,i21) 

Gross welfare this period to 
United States and rest of world 

HD, 27 6 


— 

Heading for printout. 

MN, 6 16 


"i 

Accumulator for mean values by 
period of 27 quantities listed 
in Section 1.2.7, Step Ij. 

MULT, 6 27 

MULT 

Multiplier defined in 
Section 1.2.7, Step Id. 

NO, 1 


— 

Dimension of state space next 
period, number of rows of M. 

NET, 5 


STAT(i22,... ,i27) 

Net welfare each of six categories 
this period. 

POST, 2 


STAT(il7,il8) 

Cost to producers of production 
planted this period in United 
States and rest of world. 

PHI, 2 
4 


$ 

Random term in state transformation 

PHB, 4 


— 

Expansion of 4 to maximum 
dimension of state space. 

PRC, 2 


STAT(il4,il5) 

Price of wheat this period in 
United States and rest of world. 

REV, 1 


STAT(il6) 

Revenue to United States for 
exports shipped this period. 



APL Name and Shape Algebraic Name and Description 


SIGSQ, 2 10 


Input array describing perfor- 
mance of information system. 

TOST, 1 

STAT(il9) 

Transportation cost this period. 

VAR, 6 16 

VAR 

Accumulator for variances by 
period of 27 quantities listed 
in Section 1.2.7, Step Ij. 

XB, 4 

STAT(il,... ,i4) 

State variable expanded to 
maximum dimension. 

YB, 5 

STAT(i5,... ,i9) 

Decision variable expanded 
to maximum dimension. 

The function 

SIMSTATS is direct 

implementation of the statistics 


algorithm presented in Section 1.2.7. It does not change the arrays 
MEANX and STDX as SIMGRID does, but simply prints out means and standard 
deviations of the 27 tracked quantities. 

1.3.10 Subordinate Functions 

The remaining functions are listed in Figure 1.9 through Figure 1.15. 

BUILDX builds the grid of points as described in Section 1.2.3, 

Step 3. In case the grid so built includes any points with negative coor- 
dinates, it is shifted to avoid this condition, and NOTESHIFT is called 
to print a warning. This does not happen in the normal workings of the 
algorithm, but when the statistical parameters MEANX and STDX are very 
far from their final values, it may occur. Unless it persists in itera- 
tions close to convergence, there is no problem. 

The function PIVOT is called from the quadratic programming algorithm 
MAXIMIZE, to do the pivot operation for each iteration. The left argument 
is the pivot row index, and the right argument is the pivot column index. 
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Figure 1.9 Listings of APL Functions BUILDX and CONVERGE 



n C O f-,' V C.- P: 0 - E U r n T -V 
K? Ef-COK ' verged K 

[1] 'AT ITERATION* , 4 0 

C2D (0<- ( r / I f (O f QU)-(QRREV,ORREVU) ) ) <TOLxr / I y O pTOL. 

1 : 3 ] QRREVf.a 

C 43 c?F-P:Evufau 

Z51 -iO 

V 

♦ 

V f '! o T E E H I E T [; 0 3 y 

V NOTESHIFT 

[;j 3 -jAFL.AGxNFKT 

[;2 3 NENT* 'GRID SI-IIFTED DURING R E R X C) D ' >, Q 

C33 afeagro 

V 


VR I VOTING 3 '•7 
V I RIVOT JJIRfT 

m: Fi:i3<-' 

r 2 D MATRiN[; I ; 3 ^ t< matrik[; i f 34 .matri;c[; i j J 3 
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Figure 1.10 Listing of APL Functions CONVERGEU, NOTESHIFT 
and PIVOT 
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vSTRucTURnr[;n;]'7 
^ STRUCTURE TJM4/15U 

Ci: m 44<-( \4)<>.=a4 

1 : 2.1 .i-2y<-3ri-3f i-5»i-6 

L32 2 2 f 1 0 0 1 

r-41 >>(■ 2 3 f "1 “1 0 0 1 ~1 

C51 l.7:a<-aiw 

C61 b<-ki:n 

C71 (-R 2 3 f 1 1 0 0 0 1 

C81 -iO 

[.'91 42f 1 000010 0 

Cl 01 '-.e:NR 4 5 P “1 “1 p(5fQ>9 1 0 0 0 1 0 “1 

cm <'<•■ 1 1 0 1 0 
C121 aru\u\i:21AI'< 

C131 

C 1 4 1 « C 1 ? 3 ? 3 1 C 2 ? 3 1 * 3 1 r-R-amm a [; x j 1 3 

C131 <^C1 f 55 5 K--(^amma[;i 521 

C161 r-Cl 531M:<C2531t -i:'E:LTA|:i5i3 

C171 »C1 5 5K--r'E;i.-TA[;i 5 21 

C181 -^'C4v 1531f'AC4535 11 <"GAMMA[;i5 11 

C191 <^C6 5 5y51<"(^'AMMA[:x 5 21 

C201 2 5 f 1 1 0 0 0 0 0 0 1 0 

C211 ->0 

C221 i-3;WfM44 

C2:n 4 3 f “1 "1 0 0 0 0 0 1 ~i 0 0 0 

C241 -M...7 

C251 i..5;m<“M44 
C261 ->^-S 

C271 1-6 2 4 f 1 1 0 0 0 0 1 1 

C281 2 4 f “1 “1 0 0 0 1 “1 1 

C291 LJt- 1110 
C301 '^fiJ\'-'\C21'^^^'''’ 

r 3 1 3 B u \ B :r ?■! 

C 3 2 1 « C 1 5 4 ? 4 1 ^ - f. A M M A c 1 5 2 1 

C 3 3 1 C 15 41 < ■■ 1^' C 1' 5 2 1 

C 34 1 '^^Ca 5 4 5 4 l<-f^''='MMA|;; 1 5 21 

C351 c:<- 2 4 f 1 1 0 0 0 0 1 0 

C361 -JO 

r-' 


y TABLE AU [;□!'? 

y TABLEAU 5 ZM 5 ZNM 5 XN 5 TNM 

C 1 1 ( \ ~ \ f 

[ 2 1 T h y,-l- ( \ N F: r 4. f •! C T’ ) <> , r.= X W R T' + f { C T' 

[; 3 3 tabrc ,C11(2xe:),C1 1 -hme ( np:t' , nct ) f 0 

[; 4 3 TABf TAB , IN , 1 3 ( ) y C 1 1 .> PO 

C51 tab4-tab,zn,[;i 3(6;C) ,[:il^f< 

C 61 T A B t- T A B , ( Z N M , Z N ) y [' 1 1 .T f ■! M 

C71 “>0 

y-' 


y ( 5 f ■ 0 ) y 1 


Figure 1.13 Listings of APL Functions STRUCTURE and TABLEAU 



vt. RAW [□:]'? 

V riRAW 

Ci: ap-i-t- 123456789 

C2D SITERSf-ITERS 
C3] S5IGSQ<-SIGSa 

C 43 Sr<<-( LOSS SHARE RESHARE SIGS«)*0.5 

C5D KRl 

C 6 D sr. 2 f- 4 6 fO 

172 tRLr: 1RER<-1 

CB] RRUR J 1+6tI'= *^P D 

C9] RHif. (DiMfsi.c ; IRER3 ) x ( 1 . 5 * 0 * 5 > X ( ~2+?i;’ 3: Mf 3 ) 

CIOD SD 2 Clt'iMj iRER]^sn2[: \ r>iM ^ irer]+rhi *2 

CUD 4RRLRX I 61if= E:R<-iREP;+l 
C12D -»rRLRx \ iters >.Kf.K + l 
C13D St.2<-(SD2-rITERS-0,5> *0.5 
C14D 

C15D t]P:i-<-i23456789 

C16D SDf.sr.xSD-fsr .2 

C17D K^-i 

C18D TLR} iRERf-l 

C 19D RUR ‘r'iM<-taMK[; i+ 6 tJ^F'e:p:D 

C20D RH<-rh, (r.iMtsr.[; } irer] ) x ( 1 ,5*0.5) x (“2+?i:'iMf 3) 
C21D -»Rlrx \61I'^e;r>.irep:+i 
C22D 4 fLRx \ ITERS^Kf-K + l 

V 


Figure 1-14 Listing of APL Function DRAW 



V I- s } T j V j i< ; I f Y ; Tcr j E 

nil + 

i:: 2 : t^-c k „ .TMAK)f 

C3] V,-. „ ( ^.r.IM) o 

1:4:1 

C5] i...onF- jY<..X[;;,f i;j 
Li:! Yf.Y~n: 

Z?:\ TL 5 T3<.(V/j Yc , XY) rT-„ I 
L82 -H-..00F X\ IMAX_>Tf.I + j 

L93 te^-(V/jEo , XE-) ,En, 1 

Lio: COEFFf. (HT.}., X 6 <T) + . XT+. XVAL 

LI 13 COEFFfCOEFF^TE 

<7 


r;'I...BUSLn:i'? 

V I..SUS f T f V J K J I 5 YJ TE 5 E 
L 1 :i Kf.2 J 2 ♦ 

L23 Tf.( l< , .TMAX)f If.] 

C33 Vf.j, ( 

L 4 .1 E f ( .f / L 2 1 •■•' ) -V- ^ M A 

C53 1 . 0 DP ♦ T'fML ? X] 

C63 Tf-YfE 

1.7 3 X L ; :r 3 <■- ( v/ , r « , x y > y j. 1 

LG3 -H-DOFX \ IMAK> If.I.f] 

L9.3 TEf.(V/,Eo.xE:),E,]. 

L103 COEFFf. (MATf. (E3T+. X T ) .f , XT) + , XVAL 

L 1 1 3 CDEFFUf MATf. , XVAL.U 

L 1211 rDEFFfCOEFF.fTE 

L 1 3.3 DDEFFUf-CDEFFU-TE 


Figure 1.15 Listings of APL Functions LS and LSUS 
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Function DRAW is called from SIMGRID and SIMSTATS to select the vector 
of random terms as described in Section 1.2.4 d. 

Functions CONVERGE and CONVERGED simply determine whether convergence 
has been achieved in VFS and USVFS respectively. 

Function LS calculates the least squares fit described in Step 8 
(a and b) of Section 1.2.3. Function LSUS does the same calculation for 
the United States value functions as described in Step 8 (a and b) of 
Section 1.2.5. 

Function NEWGRID calculates the new statistical parameters after a 
run of SIMGRID. It is called from ITERATE. 



2. VARIATIONS IN INPUT DATA 


As shown in the diagram of Figure 1.1, there are two categories 
of input data used for the value of information calculations. There 
are the economic parameters (discount rate, elasticities, timing 
assumptions) used in both the dynamic programming calculations and 
the simulations. Then there are the data describing the information 
systems used directly in the simulations, but also used indirectly in 
the dynamic programming calculations, since these require the statistical 
parameters (MEANX, STDX) calculated by the simulations. 

Changes in the economic parameters are easily made by replacing 
the appropriate lines of INITIAL. To correctly make changes in the 
description of either the current or the improved information system, 
one must understand how these descriptions are related to the assumed 
production estimate accuracies. 

This relationship is portrayed in Figure 3.1. The array 
describing an information system is formed by a difference operation 
from the mean squared errors in production estimates by time of 
year, together with the mean squared error of the a priori production 
estimate. Let i = 1 for the United States, i = 2 for the rest of the 
world. Let e.. , j =1, 2 , ... , 12, be the mean squared errors 

' J 

in production estimates for the United States and the rest of the 
world at two month intervals from June before planting (j = 1 ) to 
April after harvest (j = 12). We assume ~0, since regardless 

of what is published, the "market" must discover the truth as the 
supply is exhausted. 



A Priori Accuracy 



Difference 

Vector of Accuracies 
by Period of Year 

Operator 




array 




2 

Figure 2.1 Calculation of Z array— Either 
Current or Improved System 


In the case of the current system, we have obtained estimates 
2 

of the e.. 's from a study of published statistics. But these 

* J 

estimates are not based on the assumption that the true annual 

production is ever known with perfect accuracy. Built into our 
2 

's IS the assumption that a certain residual mean squared er- 
ij ( 

ror r.^ , i = l , 2 remains in the final published estimates. If 
we were to replace this assumption with the assumption of another 
value (r.' ) for the residual error, then the e. . 's would be 

I I J 

replaced with 

''ij^ ■ ^ j=l, •••. 11 ; i = l ,2, 


with 


'112 


2 


remaining 0. 




2 

For the first few periods, e.j gives the mean squared error 
Of the a priori estimate, since there is no information available 
specific to the crop about to be planted. For our wheat calcula- 
tions, this applies up to j = 6 . Now we form 

, i = 1 , 2 ; j = 1 , 10 , by 

2 _ 2 2 
°ij " ^i,j+i ~ ^i.j+2 * 

Since = e. ^ for j = l, 2, ..., 6, the first few 's are 0, 

I I i I J 

namely for j = 1 , 2, 3, 4 . 

Notice that if the residual error assumption is changed from 
r^ to r^. ' , each remains the same except for , i = 1 , 2 , 
which are simply increased by (r^') - r^. , i = l ,2 . 

In the case of an improved information system, the mean squared 
error estimates are not based on published statistics, but on 
an analysis of the sampling and measurement methods used. Thus, 
the above concept of "residual error" does not apply, except as it 
affects the a priori mean squared errors e.. for the first few j . 

’ J 

In advance of the system's measurement of the growing crop, for 
j < J (some appropriate J), we have the same e.. as for the cur- 

* J 

rent system. Therefore, if we change the residual error assumption 
of the current system from r^. to r^ ' , changes (i = l,2) , being 

increased by (r^') - r^. . The other 's are unchanged. 
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